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Abstract 



Exceptional points and double poles of the S matrix are both characterized 



> : 
o ■ 

by the coalescence of a pair of eigenvalues. In the first case, the coalescence 
^ . causes a defect of the Hilbert space. In the second case, this is not so as 

| shown in prevoius papers. Mathematically, the reason for this difference is 

the bi-orthogonality of the eigenfunctions of a non-Hermitian operator that 
is ignored in the first case. The consequences for the topological structure of 
r**~ ■ the Hilbert space are studied and compared with existing experimental data. 
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I~í ■ I. INTRODUCTION 

(N ■ 

o : 

Information on the topological structure of the Hilbert space can be obtained from a 
O-i, study of its singularities. Berry showed that geomètric phases appear when a diabòlic 
point is surrounded by varying adiabatically external parameters of a quantum system. 
Manifestations of this phase factor have been considered and proven experimentally already 
in the eighties in many different fields of physics, e.g. @. In 1994, it has been studied also 
by using microwave resonators: the sign change of the wave function has been found after a 
cyclic excursion around a diabòlic point in the space of shapes of the resonator ||. 

Other singularities are exceptional points |4j] which appear in the complex A plane of the 
eigenvalues £&(A) of the Hamiltonian H = Hq + AH\. Their positions are characteristic of 
the Hamiltonian H, once Hq and H% are given PH[7| (which both are assumed to be real 
and symmetric). The exceptional points are characterized by the coalescence of a pair of 
eigenvalues, i.e. £/(A EP ) = £ fc (A EP ). When the corresponding eigenfunctions are assumed 
to be orthogonalized in the Standard manner, it follows ipi(A EP ) = ip k (A EP ). This means, 
that ipk{A EP ) can not be normalized at A = A E p, since the orthogonality conflicts with the 
normalization requirement. As a consequence, an exceptional point is characterized by the 
fact that the rank of the associated matrix H + A E piïi drops by one at A = A E p and the 
two wave functions coalesce into one. This implies a defect of the underlying Hilbert space 

In ||, the topological structure of exceptional points is studied experimentally by using a 
microwave resonator. The exceptional point is surrounded by varying adiabatically external 
parameters of the system. As a result, the eigenvalues and eigenvectors are exchanged while 
encircling an exceptional point, but one of the eigenvectors undergoes a sign change which 
can be discerned in the field patterns. From these results, the authors draw the conclusion 
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that the exceptional points can clearly be distinguished from other topological singularities 
such as diabòlic points. 

In describing physical processes, the exceptional points lead to problems. First, the 
splitting of the Hamiltonian H into H and Hi can not be done arbitrarily. For a fixed H Q , the 
part H\ is well defined since it describes the coupling of the states of the system (described 
by Hq) via the environment (continuum of decay channels) into which it is embedded [[|. 
Secondly, the property £/(Aep) = ^(Aep) is characteristic of a double pole of the S matrix. 
Here, Sk = E^ — i/2Tk is the complex energy of the resonance state k with energy 
and width rv The S matrix describes physical processes, and no hints at all to defects of 
the Hilbert space are known at a double pole. For numerical examples see the results of 
calculations performed in a schematical model jÏ0| and for àtoms |ÏÏ|JÏ^| . 

Theoretical studies have shown that the topological structure of avoided level crossings 
is directly related to the topological structure of double poles of the S matrix being branch 
points in the complex plane |[L3|| . The transition from a double pole of the S matrix to 
an avoided level crossing by varying a parameter occurs continuously. The avoided level 
crossings are directly related to the diabòlic points 0. Thus, the topological structure of a 
double pole of the S matrix and that of a diabòlic point are related to one another. 

The problem is now the following. The double pole of the S matrix and the exceptional 
point are both characterized by the coalescence of two eigenvalues of a non-Hermitian Hamil- 
ton operator at a certain value of a parameter. Nevertheless, their topological structures 
are different: according to ||, the topological structure of an exceptional point differs from 
that of a diabòlic point while the topological structure of the double pole is related to that 
of the diabòlic point, as discussed above [jT3 |. The question arises therefore what differences 
exist between the exceptional points and the double poles of the S matrix that could cause 
their different topological structures. 

In order to find an answer to this question, the Schròdinger equation has to be solved 
in the whole function space containing everything, i.e. discrete and continuous states. By 
using a projection operator technique, an effective Hamiltonian can be derived from this 
Schròdinger equation which describes the system (Q subspace) after embedding it into the 
continuum of decay channels (P subspace) Q. Its eigenvalues and eigenfunctions are com- 
plex. The eigenvalues coincide with the poles of the S matrix. The eigenfunctions are 
related to the wave functions of the resonance states by a Lippmann-Schwinger-like relation 
Ï4J. They are bi-orthogonal. At the double pole of the S matrix, the Hilbert space has no 



defect due to the bi-orthogonality of the wave functions. The S matrix behaves smoothly 



by varying parameters also when the double pole is met QT3 . 

It is the aim of the present paper to derive the phase changes of the wave functions that 
appear after surrounding a double pole of the S matrix parametrically. In Sect. II, the 
relation between the eigenvalues of the effective Hamilton operator and the poles of the S 
matrix is scetched while in Sect. III the relation between double poles and avoided level 
crossings is discussed. The double poles of the S matrix are branch points in the complex 
plane. At these points, the wave functions of the two states are exchanged, ipi — > ±iip2- 
This causes a mixing of the wave functions in the region of avoided level crossings. In Sect. 
IV, the phase changes appearing after surrounding a diabòlic point and a double pole of the 
S matrix are derived and compared with experimental results. They agree with all data for 
isolated crossings of two states that are published in P|§. Conclusions on the topological 
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structure of the function space are drawn in the last section. 



II. EFFECTIVE HAMILTONIAN AND POLES OF THE S MATRIX 

After embedding a system into the continuum of decay channels, the discrete states of 
the system turn over in resonance states with a finite lifetime. The Hamiltonian of the 
system becomes effectively non-Hermitian with complex eigenvalues S k = E k — i/2T k , where 
the width is inverse proportional to the lifetime. 

The relation between the poles of the S matrix and the complex eigenvalues £ k can be 
derived from the Schròdinger equation 

(H-E)* E = Q (1) 

with the Hamilton operator H and the set {^e} of wave functions containing the discrete 
states of the system as well as the scattering wave functions of the environment into which 
the system is embedded. The operator H is hermitian. 
In a first step, two sets of equations have to be solved: 

_ ^cl) $ cl = (2) 

and 

E(# cc '-£KÍ (+) = o (3) 

d 

where H cl describes the system with the discrete states k and H cc ' the continuum with 
coupled decay channels c. Then, the two projection operators are defined by 

Q = E l*?X*íl P = J2[àE ití +) >tó +) l (4) 

k c J 

and H cl is identified with QHQ = Hqq and H cc ' with PHP = Hpp. The two other terms 
of H = Hqq + Hpp + Hqp + Hpq describe the coupling between the two subspaces. The 
solutions of coupled channel equations with source term 

E(# CC ' " E) (Ce (+) \ *k) = ~ {Ce (+) \Hpq\^) (5) 

d 

provide the wave functions u k that contain the coupling between the two subspaces. 

Using the completeness relation P + Q = 1, one obtains for the solution of the whole 
problem 

i N ~c 

ír| = + ' y q > . (6) 

Here, 

n k = $ fe + cj k = (i + g^Hpq) í> fc (7) 
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is the wave function of the resonance state k, Gp = P(E — Hpp) 1 P is the Green function 
in the P subspace, Cj k is determined by (H) with replaced by and 



%(e) = V2^ m\H QP \e E i+) ) = vïï (e E i+) \Hp Q \h) . (8) 

Further, is eigenfunction and £ k — E k — i/2 T k eigenvalue of the effective Hamiltonian 

W = H QQ + HqpG < p ) Hpq (9) 

that describes the system after embedding it into the continuum of decay channels. TC is 
non-hermitian, its eigenvalues and eigenvectors are complex. The eigenfunctions are bi- 
orthogonal, 

*k,i (10) 
where & k ght = <£> fc and í>). eft = p3|,|ï5| . As a consequence 

($ fc |$,) = Re«$ fe |<£ fe )); A k = (Ò k \Ò k ) > 1 

(è k \è^ k ) = i im({è k \è^ k )) = -(è^ k \è k ) ; Bf k = \{è k \è^ k )\ > o . (11) 

Using (|6|), (0) and the Lippmann-Schwinger equation for the scattering wave functions, 



one gets for the resonance part of the S matrix P,|13 



N ~c ~c' 

gires) = . y, l k l k _ (12) 

Z / 771 771 i l p v ' 

fc=l -c/ — ivjfc + 2 1 fc 

The 7^; are the coupling matrix elements of the resonance states to the continuum. de- 
scribes the resonance part of the S matrix also in the overlapping regime. The interferences 
between the resonance states are taken into account by diagonalizing the effective Hamil- 
tonian TC. Due to the unitarity of the S matrix, the j k , E k and T k are energy dependent 
functions. The relation T k = XX7fc) 2 holds only for isolated resonances. In the overlap- 
ping regime, the energy dependence of both functions is different, as a rule. For numerical 



examples see [16 



As can be seen from (|Ï2]), the poles of the S matrix are determined by the eigenvalues 
S k of the effective Hamiltonian (|J) after solving the fixed-point equations £ k — S k {E — E k ) 
|13| |. As an example, resonances of a microwave cavity are studied experimentally in the 
overlapping regime ||17|| . The results show the phenomenon of resonance trapping and are 
described well by flT2|) with the effective Hamiltonian (P). 



III. DOUBLE POLES OF THE S MATRIX AND AVOIDED LEVEL CROSSINGS 

The relation between double poles of the S matrix and avoided level crossings can be 
illustrated best by means of a simple two-level model. Let us consider the complex two-by- 
two Hamiltonian matrix 

n=( e ^ x) -^ n " ) (is) 

l u e 2 (A) - | 72 ) 
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where and 7^ (k = 1, 2) are the unperturbed energies and widths, respectively, of the two 
states. The are assumed to depend on the parameter À in such a manner that the two 
states may cross in energy at À cr when uj = 0. The two states interact only via uj which is 
assumed in the following to be independent of the parameter À (as the 7^). The eigenvalues 
of TC are 



e±--t± 

± 2 ± 



[e x + e 2 ) - - (71 + 72) 



with 



[ex ~ e 2 ) - - (71 - 72) 



+ Auj 2 



(14) 



(15) 



When F(X, uj) = at À = À cr (and uj = uj ct ), the S matrix has a double pole. 

According to (|i~5|) , F = Fr + i Fj is generally a complex number. For illustration, let us 
discuss the case with real uj. Then ei = e 2 at A = À cr and we have to differentiate between 
three cases 



F R {X,u) > 
F R (X,u) = 
Fr(X,uj) < 



'Fr = real 

= í) 
'Fr = miaj. 



(16) 
(17) 
(18) 



The first case gives the avoided level crossing in energy with an exchange of the two wave 
functions at À cr . The second case corresponds to the double pole of the S matrix. In the 
third case, the two levels cross freely in energy and the two states are not exchanged at the 
critical value A c 
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In jn|, the two cases Fr > and Fr < are studied experiment ally 
in a microwave cavity and called overcritical and subcritical coupling, respectively. The more 
complicated cases with complex u are considered in fï9|| . 

The example with real uj illustrates nicely the relation between a double pole of the S 
matrix and avoided or even free crossings of two levels in the complex plane. The double pole 
is a branch point in the complex plane. The number of these branch points is of measure 
zero, but their influence on the dynamics of quantum systems can be traced in the many 
avoided level crossings. While the wave functions of the two states are exchanged just at the 
double pole of the S matrix and are unmixed at any value of the parameter different from 
the critical one, this is not so at an avoided level crossing. In this case, the wave functions 
remain mixed in a certain range of the parameter around the critical value. This fact has a 
strong influence on the mixing of all the wave functions of a system when the level density 
is high, and different avoided level crossings appear at vàlues of the parameter inside this 
range. For the results of numerical studies see IJT3] . 

The bi-orthogonality relation ([H]) holds every where, including at the double pole of the 
S matrix. The reason is that Ak — ► 00; B l k — > 00 [Eq. (|ÏT1)1 and that (<&j£|$j) is the 
difference between two infmitely large numbers (but not their sum). This difference may be 
(for l ^ k) or 1 (for l = k). Thus, the orthogonality and normalization requirements do 
not conflict with one another and the Hilbert space has no defect at all. For the results of 



numerical studies see |L3]. 
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It should be mentioned here, that the bi-orthogonality of the {$&} follows directly from 
the non-Hermiticity of Tí. Only for the eigenfunctions of a Hermitian operator holds = 
$^ lght p_5|| . Due to the symmetry of Tí it holds $^ ft = $^ lght * for its eigenfunctions what 
results in (|Ï0|) for the bi-orthogonality relation. 

Further analytical studies |0j| have shown that the wave functions of the two states at 
the double pole of the S matrix are exchanged. It is 

(19) 

in approaching the double pole of the S matrix. This result is confirmed by numerical studies 



on làser induced continuum structures in àtoms ||12|| . 

The real and imaginary parts of the wave functions of two resonance states as a function 
of an external parameter increase limitless in approaching the double pole of the S matrix 
|Ï3fl . The sign of the imaginary part jumps at the double pole (when <£>i — > + £$2). When 
the double pole is not met by varying the external parameter, but the levels avoid crossing 
at the critical value of the parameter, the real and imaginary parts remain finite but the 
jump of the sign remains. The wave function 

<J ch = ai $1 ± i a 2 $ 2 • (20) 

changes smoothly (without any jump of the sign of its components) for a 2 — > «i at the 
double pole of the S matrix or at the critical value of the parameter where the levels avoid 
crossing. For the results of a numerical study see [Ï3[] . 



The diabòlic points are related to avoided crossings of discrete levels. They occur by 
varying two independent parameters: at the diabòlic point, two energy surfaces drawn over 
the plane of the two external parameters touch each other at one point forming a double 
cone. 



IV. GEOMÈTRIC PHASES 

Let us now consider the geomètric phases appearing after encircling a diabòlic point and 
a branch point in the complex plane (double pole of the S matrix), respectively. In any 
case, the paths of encircling are characterized by the value F, Eq. flUp, which vanishes only 
at the branch point in the complex plane. Most interesting are states whose eigenvalues are 
near to the real axis. We can restrict our discussion therefore to real uj (see Sect. III). 

For encircling the diabòlic point or the branch point in the complex plane, two external 
parameters have to be varied. In the experiment the diabòlic point is surrounded by 
varying the shape of the microwave resonator by means of two parameters but leaving the 
coupling strength to the antenna unchanged. Since the two levels considered avoid crossing, 
the whole path of encircling the diabòlic point is in the overcritical regime. That means, 
the critical value of the parameter is passed twice, on the way forth as well as back, under 
overcritical conditions, and the wave functions are exchanged each time when the critical 
value of the parameter is reached. This is not so in the experiment || where one of the 
two parameters is the coupling strength of the cavity to another one. Therefore, the critical 
value of the parameter is passed on the path of encircling the exceptional point (or branch 
point in the complex plane) only once under overcritical conditions. The other part of the 
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path is in the subcritical regime where the wave functions are nowhere exchanged (see Sect. 
III). 

In detail: The diabòlic point is surrounded in the experiment || in the regime of overcrit- 
ical coupling along the whole way of encircling and À cr is passed twice in opposite directions, 
(i) $ fc -> -ièf, -> +zí> fc , i.e. 

{$1, í> 2 } -> {-zí> 2 , +zí>i} (21) 

and (ii), on the way back, — > — i <3>fc; — > +i $z, i.e. 

{-ï$ 2 , +zí>i} -> {-l>i, -<§ 2 } . (22) 

The phase change occuring after one surrounding the diabòlic point is therefore 

{$i, $2} -> {-$1, -$2} • (23) 

This corresponds to the geomètric phase discussed by Berry JT[j. 

The way of encircling the branch point in the complex plane itself passes from a region 
with overcritical coupling at À cr to another one with subcritical coupling at A cr . An exchange 
of the wave functions takes place only at overcritical coupling where the resonances avoid 
crossing. Thus, a first full surrounding gives 

{$!, <l 2 } -> {-iè 2 , +1^} (24) 

and a second one (in the same direction) 

{-iè 2 , +2$!} -> { + $!, +<J 2 }. (25) 

That means, surrounding the branch point in the complex plane twice restores the wave 
functions &k including their phases. This corresponds to the result obtained for surrounding 
the diabòlic point twice. In both cases, the wave functions including their phases are restored 
after a second encircling in the same direction: 

{$!, í> 2 } => {<£>!, í> 2 } . (26) 

Encircling the branch point in the complex plane in the opposite direction gives 

{<£>!, í> 2 } -> {+1Ò2, -z$i}. (27) 

Since the experiment [U is not sensitive to the possible occurence of a phase i of the wave 
function, the results ( p4"| ) for one loop with a certain orientation of the path and fl2"7j) with 
the opposite orientation of the path agree with the experimental data given in ||. There 
are no experimental data in || for the phase changes after a second loop. 

An experimental study of interferences between atòmic levels in a làser field is expected 
|T9| to allow conclusions on the phase changes, including those after a second loop. 
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V. CONCLUDING REMARKS 



In the present paper, the phase changes occuring after encircling parametrically an iso- 
lated diabòlic point and an double pole of the S matrix (branch point in the complex plane) 
are calculated. The results are shown to agree with all experimental data that are published 
in 

The results of @ point to the interesting fact that the phase changes after surround- 
ing higher-order degenerades are more complicated than those obtained after encircling a 
diabòlic point. This result has given rise to further theoretical studies, e.g. [fJÜ] . 



The experimental results |J are interpreted by the authors on the basis of exceptional 
points. This interpretation leads to the conclusion that an exceptional point can clearly 
be distinguished from other topological singularities such as diabòlic points. The authors 
claim the following: encircling the exceptional point a second time completely with the same 
orientation one obtains {— — while the next complete loop yields {— and only 
the fourth loop restores fully the original pair The authors show experimental 

results only for one complete loop. No data are given for two or more loops. 

The appearance of a phase change of both wave functions after a second loop around 
the exceptional point, <£>fc =>- —&k', &i =>• suggested in ||, does not agree with the 

result ( P^ ) obtained for a second complete loop around a branch point in the complex plane. 
According to (|26|), the original pair is restored already after a second complete loop 

when it is completed with the same orientation. This result coincides with that obtained for 
a second loop around a diabòlic point. It is an expression for the fact that diabòlic points 
and branch points in the complex plane are related to one another as discussed in this paper. 

The results for one loop can not differentiate between the two interpretations since the 
experiment is not sensitive to the possible occurence of a phase i in the wave function. 
It can therefore not be concluded from the published experimental data whether or not 
the topological structure studied in || is different from that of a diabòlic point. Further 



experimental studies are necessary, maybe on àtoms in a làser field as suggested in f[Ü 
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